function f = plotf2(p0, gamma, lambdaLower, lambdaUpper)

m = 500;
h = (lambdaUpper- lambdaLower)/m;
for i = 1:m
	lambda(i) = (i-1)*h + lambdaLower;
    f(i) = tanh((lambda(i)^.5 - (lambda(i) + 2*p0*gamma)^.5)/gamma) - ...
         lambda(i)/(9.8+gamma*(lambda(i)^.5));
    df(i) = (sech((sqrt(lambda(i))-sqrt(lambda(i) +2*p0*gamma))/gamma))^2 * (1/sqrt(lambda(i)) - 1/sqrt(lambda(i) + 2*p0*gamma))/(2*gamma) - ... 
(9.8 + .5*sqrt(lambda(i))*gamma)/((9.8 + gamma*sqrt(lambda(i)))^2);
end

plot(lambda,f);
figure;
plot(lambda, df);